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Abstract: Cosmic structures at small non-linear scales k > L ~ 0.2h Mpc -1 have an 
impact on the longer (quasi-)linear wavelengths with k < L via non-linear UV-IR mode 
coupling. We evaluate this effect for a ACDM universe applying the effective fluid method 
of Baumann, Nicolis, Senatore and Zaldarriaga [1]. For k < L the ACDM growth function 
for the density contrast is found to receive a scale dependent correction and an effective 
anisotropic stress sources a shift between the two gravitational potentials, setting (j) - ip 7^ 0. 
Since such a situation is generically considered as a signature of modified gravity and/or 
dark energy, these effects should be taken into account before any conclusions on the dark 
sector are drawn from the interpretation of future observations. 



1. Introduction 



The standard ACDM model of the universe is characterized by a specific evolution of 
the density contrast 5p/p = 5 oc aD(a) i.e by a specific growth function D{a) for the 
matter perturbations, and by the equality of the two gravitational potentials 4> and ip as 
is implied by the lack of anisotropic stress in the CDM. In contrast, evolving dark energy, 
or a modification of Einstein gravity, would generically give rise to an anisotropic stress 
as well as to modifications in the Poisson equation [2], which would manifest themselves 
in a different functional form for D{a). Probing the growth function for different values of 
scale factors a would thus be invaluable to determine the dark energy related cosmological 
parameters. It is for this reason that the three-dimensional mapping of the structure of 
the universe by future weak lensing surveys such as EUCLID [3] is very much expected to 
significantly constrain dark energy and modified gravity models. 

At a practical level the determination of the growth function proceeds from the as- 
sumption that the universe can be considered as a collection of ideal fluids with some small 
perturbations. In this setting, an accurate (scale independent) parametrization for D reads 



where fi m (a) = HfamCT^ / H (a) 2 and 7 (= for ACDM) depends e.g. on the equation 
of state parameter of the dark energy. 1 However, structure formation itself feeds back on 
the stress-energy and hence on the effective behavior of the cosmological fluids. As has 
been pointed out in [1], this effect can be studied by integrating out the short-wavelength 
perturbations, obtaining an effective theory for the long-wavelength universe that has an 
equation of state different from the homogeneous background and, moreover, is no longer 
strictly ideal but is also characterized by a viscosity parameter. Very roughtly, one can 
think that this effect is caused by the motion of small-scale lumps of matter and the tidal 
effects of longer wavelength perturbations on this motion. However, one should also note 
that the scales that have virialized can actually be shown to decouple completely from the 
large-scale dynamics at all orders in the post-Newtonian expansion, apart of course from 
contributing a correction to the energy density [1]; in effect, for a distant observer virialized 
systems behave as if they were point particles. 

In the present paper we follow this approach and study the effective cosmological fluid 
of the long-wavelength perturbations in a pure ACDM universe. We point out that the 
naive expectations for e.g. the growth function in the linear regime are violated at some 
level due to the effective pressure and anisotropic stress generated from the non-linear mode 
coupling between short and long wavelengths. This effect is generic and scale-dependent 
and should therefore be taken into account before any conclusions about the properties of 
dark energy are drawn from the interpretation of observations. On the other hand, it could 
also be used as a consistency check of the concordance model. 

x It is interesting to note that scale dependent growth functions are physically more realistic for dark 
energy models [5]. 



[4] 




(1.1) 
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We proceed with our investigation as follows: in the next section we briefly review 
the ingredients of the effective fluid approach to the non-linear coupling of long to short 
scales [1] that we need to derive our results. Using the methods described in section 
2, we compute in section 3 the evolution of perturbations in this effective fluid and the 
corresponding corrections to the linear ACDM growth function, as well as the so-called 
gravitational slip: the divergence of the two gravitational potentials (f> — ip. In the process 
we also provide an exact analytic expression for the linear ACDM growth function which, to 
our knowledge, has not appeared in the literature previously. We summarize and conclude 
in section 4. 

2. The Effective fluid 

As matter perturbations in the universe grow under the influence of gravity there finally 
comes a point for a given scale where the density contrasts become larger than unity. At 
this point naive perturbation theory breaks down and other techniques need to be used 
in order to follow the further evolution of perturbations, eventually having to resort to 
numerical simulations. In our universe scales with k > k n \ ~ 0.2h Mpc -1 have entered the 
non-linear regime. 

Apart from the difficulty of following the evolution of non-linear perturbations, another 
issue is the inevitable coupling between different scales which is absent in the linear regime. 
In particular long and short scales influence each other and it is expected that the formation 
of non-linear perturbations might have an impact on larger, linear scales through UV-IR 
coupling. It has even been suggested that even the observed acceleration might be due 
to such non-linear effects, and this possibility has recently spurred a lively debate in the 
literature [6] . Regardless of whether such an explanation for the observed acceleration is 
feasible, the UV-IR coupling is always present and should be quantified, especially given 
the accuracy of forthcoming observations which promise to map the universe and its history 
with unprecedented accuracy. 

The authors of [1] address this non-linear mode coupling in a way that is simple 
and physically intuitive. By smoothing out short scale non-linearities they are led to 
an effective approach for longer wavelengths where the CDM fluid on scales k < k n \ is 
replaced by an effective fluid which, unlike the underlying CDM matter content, possesses 
anisotropic stress and pressure. In more detail, the basis of this approach is an effective 
energy momentum tensor, including the effects of gravity, which is smoothed out on a scale 
L < k n \. To second order in the derivatives of the gravitational potential and peculiar 
velocities, this procedure adds to the rhs of the Einstein equations the terms 



The subscripts / and s refer to "long" modes with k < L and "short" modes with k > L 



[r%] L = -[pvy s ] L - 



[t 1 3 ]l = [pvy s ] L - 



_ [di<t> s di<l>s]L ~ jW?j>s\L 

8vrGa 2 

[dicpsdicps]^ - 2[di4> s dj(j) s ] L 
8vrGa 2 



+ 0[(dv l f/L\{d<p l ) 2 /L 2 ] (2.1) 



+ 0[{dv l f/L\{d ( j H ) 2 /L 2 } . (2.2) 
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respectively. By the notation [X]l we mean 

[X] L (x) = J d 3 x'^ L (|x - x'|)X(x') , (2.3) 

where Wl(|x — x'|) is a window function that eliminates modes with k > L. Note that 
since we are only considering terms linear in metric perturbations (but not their spatial 
derivatives) and quadratic in velocities, the expressions in t^ u coincide to this order with 
proper volume averages. For k <C L the third terms on the rhs of (2.1) and (2.2) are 
suppressed. Note also that the expressions for (2.1) and (2.2) contain short wavelength 
modes with k > L but the Fourier modes of [t°o]l and \t % j\h themselves are only defined 
for k < L. For the discussion that follows it is the spatial part [t 1 j]l that will be most 
important. 

According to [1], the effective stress tensor can be expressed as 

[tS-]l(x) = <[t*,-]l> + Ar^(x) + aS-(x) • (2.4) 

with the first and second terms admitting an effective description to lowest order in 
{k/kNh) 2 as the stress tensor of an imperfect fluid with pressure and anisotropic stress 

{[r i j ] L ) + Ar i j = {V + AV)5) + ^ j , (2.5) 

with S l j = 0. In particular, the pressure terms can be written as 

V = l([r i i ] L ) = wp, (2.6) 

AV = c 2 sP 5 (2.7) 

where w and c 2 s are the equation of state parameter and the sound speed squared, while 
the scalar anisotropic stress is given as 

^ = V - #v , (2-8) 

where i] > is the coefficient of shear viscosity. Note that 

<[tS-]l> = , (2.9) 

where f l j = T l j — J*/3r'/, so there is no zero mode for the anisotropic stress. It is 
important to stress at this point that although the universe only contains CDM, the effective 
description on long wavelengths involves both pressure and shear viscosity. It is customary 
to define the anisotropy scalar a as 

' = -^-1^ = 4.^ <™> 

where we have defined the dimensionless parameter c\ is = 1rfH.j?>p. 

The possibility of expressing the effective stress-energy tensor in terms of long wave- 
length perturbations such as in eqs (2.7) and (2.8) is a manifestation of non-linear mode 
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coupling relating the short wavelength modes in the definition of [t^Jl to modes with 
k < L. The effective fluid is then described by two parameters c 2 s and c 2 is which can be 
determined by using the definition (2.2) via 

Cms ~ p k 2 <V-u,V-u,>- 

The tidal effects of long wavelength gravitational perturbations on shorter wavelength 
modes produce correlations of the long wavelength 5 and v to the short wavelength modes 
in [T l j]i, making c 2 and c 2 is different from zero. 

Finally, there will also be a stochastic part, given by the a l j term. It encodes deviations 
from the average value that are uncorrelated to long wavelength variables and acts as an 
external source. The variance of these terms can be obtained via 

<«>M = ([r'M^ik) ~ ([t^W^Wl) ~ (Ar^AA). (2.13) 

Summarizing, we see that the UV-IR coupling can be described by the use of an 
effective theory where the short scale fluctuations act as an effective imperfect fluid on 
longer wavelengths. The Einstein linearized equations with Ar l j + a % j added to the rhs 
can then be used to follow the linear response of long wavelength perturbations to short 
wavelength non-linearities. The effective fluid is characterized by two parameters: c 2 and 
c 2 is . The correlators on the rhs of (2.11) and (2.12) can be determined from small-scale 
N-body simulations using the definitions (2.1) and (2.2) for [t°o]l and [r l j]i. Once these 
parameters are determined, no further reference to the short scale dynamics is needed. 
Alternatively, the parameters of the effective fluid could be considered as free parameters 
to be fitted from observations - see section 3. 



3. Perturbations of the effective fluid in ACDM 
3.1 Density contrast and gravitational slip 

We can now proceed to evaluate the effect of short wavelength fluctuations on longer 
wavelengths using the approach outlined above. We start with the Einstein equations 
relating the gravitational potentials with the components of the effective energy momentum 
tensor. The 00 and Oi equations give 

k 2 4> + mU + %A = -^o 2 — s (3.i) 

^ + n ^ = ^ H ^(l + W )v, (3.2) 
while the traceless part of the ij equation gives 

k 2 (4>-^ = ^H 2 ^(l + w)a. (3.3) 
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We can further use the energy-momentum, conservation law V ^T^ 11 = which results in 

5 = {l + w) (k 2 v + 30) - m - w5j , (3.4) 

v = - (l-3w + — —\uv-—^— — + a-^ (3.5) 

V l + wj l + Wp m 

with v the peculiar velocity potential: Ui = ih L v. 

According to the discussion in the previous section we can express the pressure per- 
turbation and the anisotropic stress of the effective fluid in terms of long wavelength per- 
turbations 6 and v and a stochastic part uncorrelated with them 

#P = c 2 s p m 5 + c 2 s p m ai (3.6) 

a = c v is ~}r~ + c ^ s " 2 ~ c ^ is6 + c2visOL2 ' ^ 3 ' 7 ' ) 

where we have used the linear perturbation theory relation 5 ~ —k 2 v/Tl for the (long 
wavelength) velocity. The quantities a\ and a 2 encode the (dimensionless) stochastic 
fluctuations of 5P and a. We stress that we ignore here the non-linearities of the long 
wavelength perturbations 2 . 

We can now obtain an equation for the evolution of the long wavelength density con- 
trast 5. To simplify the result we make the following assumptions: 

1. We are interested in following perturbations on scales smaller than the horizon so 
that k 2 >H 2 . 

2. We will also assume that c 2 and c 2 is are small enough such that c 2 ^ < 1. This 
requires c 2 < 1CT 5 {^f^j where L is the smoothing scale we use to define the long 
wavelength sector and k D \ ~ 0.2 h Mpc" 1 is the scale of non-linearity. 

3. The time scales of evolution on scales k > L is assumed not to be much faster than 
cosmological time scales, ie J^(ca) ~ Hca and that 



Assumptions 2), 3) and 4) can of course be checked once the parameters have been calcu- 
lated from first principles. 

Given the above we obtain for the density perturbation to leading order in c 2 

S + nS- ^tf 2 — 8 = -k 2 c 2 5-k 2 (c 2 ai - c 2 vis a 2 ) , (3.8) 

which reduces to the standard equation when short-scale non-linearities are ignored ie 
c 2 — > and c 2 is — s> 0. It is useful to express the above in terms of the scale factor a. Using 

ij = an l we have for 6 ( a ) 



d\^ ( d , v , 2 d 3fH \ 2 n m 1 fc 2 p . 2 2 v 



(3.9) 



2 We could have included terms quadratic in the gradients of the potentials and the velocities, studying 
non-linearities on long wavelengths via second order perturbation theory 
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We denote the solution in the absence of mode coupling (rhs is taken to be zero) by 5 

^ + t k7{+ atf-2Uj ^ = °- (3 - 10) 

Then, to leading order in c 2 , the solution of the long wavelength density contrast in the 
presence of small scale non-linearities is 

oo 

dx G(a, x )-^ ^ x y ( c ^(» + c s"i ~ c lis a 2) , ( 3 - n ) 
where G(a, x) is the appropriate Green function satisfying 

*&a. ,) + l± ta« + ^G(«, x) - | (f ) 2 ^G(«, x) = ,„<• - ,) . (3.12) 

The solution to (3.10) that grows like 6 oc a at early times (during matter domination) 
reads 

a 

5{a) ~ Wn2 H ^ m — J Uixf' ( } 



where aj ra is the scale factor at which we start the computation and <5j n is the density 
contrast at that time. In a ACDM universe we have 

m \ 1/2 

H = H a(-^ + tt A \ , (3.14) 
and the integral can be performed to give 

a 

du 2 



3/2 n 3/2 



where we have defined F 



a 5/2 F{a), (3.15) 



with 2-^1 the hypergeometric function. Thus, 



5(a) = + ^aA V2 aF(a), (3.17) 



and the ACDM growth function for the density contrast is 



D(a) = ( 1 + ^a 3 F(a) . (3.18) 



rn 



The Green function for the problem (satisfying homogeneous initial conditions at a = 0) 
is found to be 

G(a, x) = Q(a - x) 2 —- x 2 U{x) 2 ^^- { a 5/2 F(a) - x 5/2 F{x)) . (3.19) 
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Figure 1: ACDM linear growth for density contrast. Dashed line is for the Einstein-de Sitter 
Universe: 5 ex. a. 



Using the above, we can now obtain for the solution for the density contrast from 
(3-11) 

= ( 1 ~ K$ E i( a )) " §5(a, a) , (3.20) 



where 



*w-i/*( i+ ^r^( i -er^)^ ^ 



S(a.n) - -— ( rFui) ( l + ^a 3 V /2 / dx ( 1 - (-) *'* , (3.22) 



and 



5Q m \ £l m J J \ \aJ F(a) J 6 

o 

with a = a.\ t-ol2 and c 2 s the speed of sound today. For the difference of the two 

gravitational potentials we get 

n m [l-l (l + ^a 3 ) 1/2 F(a) 
- V) = f H%cl ls —^- J ^(a) + |^o 2 — cL«2 • (3.23) 



aF(a) 



3.2 The matter powerspectrum and weak lensing 

In the previous subsection we saw how the impact of UV-IR coupling on the growth of 
perturbation can be described in terms of an effective quasi-linear theory on scales larger 
than L. The main results where equations (3.20) for the density contrast and (3.23) for 
the gravitational slip. From (3.20) we immediately obtain 

/ h 2 2 \ _ _ h 2 2 2 2 

(6 k 5 p ) = 1 - 2-^E l (a) (5 k 5 p ) - -^^(S k S p ) , (3.24) 
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were we have used — k 2 (j) ~ ^-^HqS. To simplify our formulae we will make the assumption 
that 

(S k S p ) = (2vr) 3 «5(k + p) P(a, k) P s (k) , (3.25) 

where Ps{k) is the density contrast power-spectrum: (S^Sp) = (27r) 3 <5(k + p) P$(k). In 
particular this entails that different modes of the stochastic fields are uncorrelated for 
different momentum magnitudes, at least to this order in c 2 . Given this assumption we 
can write 

<W = (l - 2 ^pi(«) - > (3-26) 

where 

P(a,k)= { *^ (3.27) 
denotes the ratio of the stochastic to the density power-spectra. Using 



^Ka-k) * I (([ri] 2 Lk ) - (N] Lk ) 2 ) (3.28) 



we obtain 



where 



4P - E 2(«);^72 (3-39) 



^-^'^('-©"'IDt' (3 ' 30) 



.3 / dq H*P S (qf 



^ L4xl ° Jk^^T- (3 - 31) 

and T(k) is the transfer function for ACDM (see eg. [7] eq. 6.5.12). To obtain (3.29) we 
used the fact that since k is in the linear regime P\m(k) = 2tt 2 5 2 H ^T{k) 2 j^fp an d Sh is 
the perturbation amplitude at the current horizon scale: 5h = 4.9 x 1CT 5 [8]. 

We have now reached our main conclusion. We see that we obtain a scale dependent 
correction to the pure ACDM growth function 

D{a)^Q{k,af/ 2 D{a)^ (^1 - 2^ cjE^a) - ^^^^^^ia)^ D(a) . (3.32) 

This may have interesting ramifications for weak lensing observations, which probe the 
potential kikjpip m (((f> + V0k(<A + VOp) [9]- Because of the UV-IR coupling, this will be also 
modified from its pure ACDM form. Defining ip = ncj) + A we obtain from (3.23) 

15 ^ 2 i-i(i + ^ 3 ) 1/2 *>) 

(i+fe- 3 ) 372 ^) 

and 



V — c vis , x 3/2 (3.33) 



A = ^ffcLa 2 (fc). (3.34) 
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The weak lensing potential can then be calculated to be 

kikjpipm^ + V>) k (0 + V0p> = [4(1 + |)Q(fc, a) + ctMk) 

2 



/3 2 f2 m \ fcjfcj pip m ~ ~ . 



with ^(k) = (ct2kQ ; 2-k}/(^k;^-k} ~ l2{k)- We can ignore terms that are not enhanced by 
factors of k, obtaining 



M:;PiPm<(0 + VOk(0 + ^) P > ^4Q(fc,a) (^o^)'^^^) • (3.36) 

In our case the weak lensing potential directly probes the modified growth function with 
contributions from the effective anisotropic stress being subdominant. 

An accurate determination of the magnitude of the parameters c 2 and 72 controlling 
the correction Q — 1 is beyond the scope of this paper since it would require calculations 
in the non- linear regime in two distinct ways: 1) Determination of 72 requires knowledge 
of the non-linear powerspectrum. 2) The effective sound speed (? s is a distinctly non-linear 
phenomenon which at lowest order depends on the three-mode coupling between one long 
and two short-wavelength modes - see (2.11) which is zero in linear theory. Therefore, it is 
more straightforward to consider that a scale-dependent correction to the ACDM growth 
function due to non-linear mode coupling would have the form 

W) _( 1 _ 2 ^ W __^^ M ), (3.37, 

where fi\ and ^2 are parameters to be fitted by observation. 

However, before closing this section let us make the following observations. We expect 
that as non-linearities grow, so will the sound speed c 2 . We can thus obtain an upper bound 
on the value of the functions El(a) and E2(a) today (a = 1) by setting c 2 s {a)/c 2 s — > 1. This 
gives 

El(a = 1) < 0.26 , E2(a = 1) < 0.1 (3.38) 

Furthermore, simply extrapolating the linear theory powerspectrum into the non-linear 
regime to calculate 72 and using 2nd order perturbation theory [10] to evaluate c 2 s we 
obtain 

c s~%tj! / dK K T 2 ( K ) -0(1) x4.2 x KT 6 , (3.39) 

2 H Q J 

V2k nl /k C q 

72 ~ $h J dKT 4 (n) ~ 10" 13 (3.40) 

We thus see that at the non-linear scale k n \ ~ 0.2 h Mpc -1 linear theory extrapolations 
predict corrections to the growth function of the order of 10% and 1% respectively for the 
two correction terms in Q (3.32), with the second term growing faster than the first with k. 
We stress again that these estimates are based on linear and second order extrapolations 
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in the highly non-linear regime. However, they could be considered indicative of a possibly 
measurable effect towards the end of the linear regime, accessible through forthcoming 
combinations of weak lensing and deep redshift surveys such as EUCLID [3] . 

4. Summary and Discussion 

We have considered the influence of short wavelength matter fluctuations on longer wave- 
lengths through non-linear mode coupling in ACDM, applying the method of [1] in this 
case. Integrating out the short scales gives rise to an effective long-wavelength fluid cou- 
pled to gravity that has slightly different properties as compared to the pure ACDM. The 
effective fluid is viscous and exhibits a small pressure as well as an anisotropic stress, all of 
which are in part correlated with the longer wavelength perturbations and in part purely 
stochastic. 

We have discussed the nature of these corrections and derived their contributions to 
the growth function and the weak lensing potential of standard ACDM cosmology. In the 
process we also provided an analytic expression for the growth function of pure ACDM 
which, to our knowledge has not appeared in the literature before, equation (3.18). The 
contribution from the induced anisotropic stress turns out to be negligible for present obser- 
vational purposes, but the modification of the growth function Q, which is scale-dependent, 
could be of importance. It contains a term that scales as k 2 and is due to the perturbation 
of the effective pressure, while a k 3 /T(k) 2 term appears because of the stochastic nature of 
small scale fluctuations that are uncorrelated to long wavelength variables and effectively 
act as an external source in the evolution equation of the long-wavelength perturbation. 
The corrections are characterized by two parameters, c 2 and 72, which can in principle 
be calculated but, perhaps more realistically, can also be fitted from observations. We 
extrapolated perturbation theory to obtain estimates for these parameters which should 
more appropriately be considered as lower limits; a more precise estimation would require 
one to go beyond perturbation theory into the non-linear regime of structure formation. 

We would like to end by noting that searches for modifications of gravity and/or 
dynamical dark energy focus on the the modifications of the growth function compared to 
its ACDM form as well as a non-zero difference in the two gravitational potentials. However, 
as we have discussed here, such modifications arise at some level even in standard ACDM 
because of the non-linearity of gravity which couples different scales. In this case, the 
difference in the two potentials seems too small to be observationally relevant but scale 
dependent corrections to the growth function might be detectable. These effects are always 
present and should be taken into account in future surveys before any conclusions are drawn 
with pertaining to the dark sector. Furthermore, such effects could be used as consistency 
checks of the standard cosmological model, provided that a more accurate determination 
of the parameters of the effective fluid is achieved. 
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